Variable range hopping and quantum creep in one dimension 
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We study the quantum non Unear response to an appUed electric field £ of a one dimensional 
pinned charge density wave or Luttinger liquid in presence of disorder. From an explicit construction 
of low lying metastable states and of bounce instanton solutions between them, we demonstrate 
quantum creep v — e~'^^^ as well as a sharp crossover &t E = E* towards a linear response 
form consistent with variable range hopping arguments, but dependent only on electronic degrees 
of freedom. 

PACS numbers: 



I 

CZ3 



I 

o 
o 



> 
m 
cn 

m 
o 
m 
o 



I 

o 
o 



X 



Computing the response of a disordered elastic system 
to an external driving force is a long standing problem. 
This is of theoretical importance and also relevant for a 
host of experimental systems, both classical and quan- 
tum. For classical systems, typical experimental realiza- 
tions are domain walls^Q, |3| and vortex lattice in type 
II superconductors 0, 0|7 Pinned quantum crystals are 
charge or spin density waves 0, Wigner crystal in two 
dimensional electron gas Q and disordered Luttinger 
liquids In the absence of quantum or thermal fluc- 
tuations disorder leads to pinning or localization. It was 
initially believed that thermal activation over barriers be- 
tween pinned states would result 9] in v{F) ~ aF albeit 
with an exponentially small mobility a. However, the 
glassy nature of such disordered elastic systems leads in- 
stead to divergent barriers and to a non linear response 
[lollTll of the form v = exp(— known as creep 

In quantum disordered systems barriers between the 
many metastable states can be overcome by thermal and 
quantum activation. Determination of the relation v{F) 
is thus an even more difficult and mostly open question. 
Two main issues arise: (i) does one recovers a quantum 
creep formula at T = when the system can unpin via 
quantum tunnelling over barriers; (ii) does one recovers 
linear response at T > 0, v{F)/F — > a and what is the T 
dependence of the conductivity a. Although these ques- 
tions have been answered in details via controlled instan- 
ton calculations for pure systems such as the Sine-Gordon 
model 0, 0, Q , with and without dissipation, no con- 
trolled method has been found for the disordered prob- 
lem. Results were obtained using physical arguments for 
very disordered electronic systems [l^ • The renormaliza- 
tion method used for creep in classical systems was 
extended to quantum problems 0| , but suffers from the 
same limitations The conductivity of charge den- 

sity waves was studied by Larkin and Lee |l9l| . but only 
in a strong pinning regime considering tunnelling around 
single impurities. 

In this paper we study the driven quantum dynamics of 
a pinned ID charge density wave or of ID interacting elec- 



trons (Luttinger liquid (LL)) in the localized phase, per- 
forming a controlled calculation of the tunnelling rates. 
It is known that this system renormalizes to strong dis- 
order where the (classical) ground state can be found ex- 
actly and low lying kink-like excited states constructed. 
We then study instanton (bounce) solutions and estimate 
the semiclassical tunnelling rate between these states, in 
presence of an applied (electric) field. This demonstrates 
a quantum creep law v = e~'^^^ ' at zero temperature. 
At small non zero temperature we show that a sharp 
crossover occurs between quantum creep for E > E* 
and linear response for E < E* . The temperature de- 
pendence of the conductivity is of the form a oc e~'^^'^ ^ 
consistent with Mott's variable range hopping (VRH) ar- 
guments (2^. Applied to the Luttinger liquid, this ex- 
tends in a precise way the validity of VRH formula to 
interacting electrons in d = 1. Note that here contrarily 
to standard VRH arguments the prefactor c of the tem- 
perature dependence in the exponential is determined by 
the electronic degrees of freedom, and is not dependent 
in an essential way on coupling to other degrees of free- 
dom such as phonons. This leads to quite different energy 
scales for c than the standard VRH mechanism. 

We consider the Hamiltonian of a charge density wave 
where the density has a sinusoidal modulation 



where 
phase I 



p{x) = po cos{Qx - (/)(x)) (1) 

i is the phase of the charge density wave. The 
obeys the standard phase action 
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where u is the velocity, y = ut and /? the inverse of 
the temperature. Furthermore the system has a short 
distance cutoff (lattice spacing) a. The disorder is mod- 
elled by a random potential V{x) coupled to the density 
by i7 = — / dxV{x)p{x). Assuming that (j) varies slowly 
at the scale we can only retain the Fourier compo- 
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nents of V{x) close to Q |22|- This leads to the action 
S,i./h = -\j ^^^j)e'W-.«)-CW)+h.c. (3) 

where we represent the disorder with a random amplitude 
A and phase C, which are both slowly varying variables. 
For a Gaussian disorder initially, the disorder ^(a;) — 
A{x)e'^^^^'' obeys £,{x)i{x')* = D5{x — x') other averages 
are zero. Adding an external electric field E to the system 
adds to the action 



(4) 



Se/^= / dxdyE(f>{x,y) 



with E = Epo/{Quh). The action 12I4II also describes 
a LL in presence of disorder l22j|. In that case 
Q — 27r/9o = 2kF where kp is the Fermi wavevector for 
fermions. K is the standard Luttinger parameter that 
describe the interactions effects {K — 1 for noninteract- 
ing electrons and K < 1 for repulsive interactions). Our 
study thus directly gives the conductivity of disordered 
LLs. In that case the pinning of the phase variable 
corresponds to the Anderson localization of the system. 

At T = the disorder is a relevant variable. It pins the 
phase (j). In the ground state the phase (f> varies by a quan- 
tity of order 2tt over a distance ^]nr which is the pinning 
length of the charge density wave [2l| or the localization 
length in the presence of interactions for the interact- 
ing particles 112 • To determine the dynamics of this 
model, we renormalize the system up to a point where the 
disorder is of order one. Since we are interested in the 
limit of very low temperatures and fields, we can renor- 
malize the action in the absence of E and at T = 0. The 
flow in that case is well known and we do not re- 

produce it here. The disorder D scales to strong coupling, 
and the parameter K decreases. We stop the flow at the 
lengthscale I* for which A ~ 1. At that lengthscale the 
disorder being of order one, the pinning length is of the 
order of the lattice spacing a. The original localization 
length of the system is thus fioc = ae' The electric held 
is also renormalized and becomes -^r^ = E{1*) = i?e^' 
and time and space are rescaled by a factor e^' . In what 
follows we denote with a star the renormalized quantities 
at the scale ^loc. Since u also renormalizes one can ab- 
sorbs this renormalization by rescaling the time by u* /u 
which changes u — *■ u* in all the above expressions. 

To study the dynamics we consider H2l4f) with the 
renormalized parameters. Although we stopped the flow 
when the disorder is of order one we assume that we are 
truly at strong disorder and can thus consider that the 
amplitude of the disorder is very large. The main effects 
thus come from the fluctuations of the random phase of 
the disorder. In order to perform a semiclassical approx- 
imation for the dynamics one must first determine the 
(classical) ground state of the renormalized system for 
E = 0. The disorder being time independent the action 



is minimized by dr<j> = 0. It is convenient to go back to 
a lattice description. The energy on the lattice is 
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u*h ~ 2TrK*a ^ 
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[(0,+i-,^,)2-A*cos(0,-C^)] (5) 



with N = L/^\oc and we take — tt < Q < tt. Since the 
renormalized disorder A* in (jSjl can be considered to be 
large (A* ~ 1), to minimize the cosine term one needs to 
take '/>i = Ci + 27rni where rii are integers. The energy 
becomes 
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with /, = (0+1 - G)/(27r), 1- Contrarily to 

higher dimension, here in d = 1 there is no frustration 
and one can minimize the action for all bonds simultane- 
ously (i.e. all pairs Arii — ti^+i — n^) by choosing [j^ : 



mo 



(7) 



where mo is an arbitrary integer and [x] denotes the clos- 
est integer to x. [x] = if —1/2 < a; < 1/2 and [x] = 1 
(resp. [x] = -1) for a; > 1/2 (resp x < -1/2). The 
values [fi] thus completely characterize the ground state. 
Here one takes the Q uniformly distributed, hence the 
perform a random walk and the ground state has rough- 
ness exponent 1/2 (i.e. 0(i) scales as (j}{L) ^ L^/"^) in 
agreement with other calculations [isj . 

In presence of the electric field E any one of these 
ground states (with toq fixed) become metastable since 
the phase wants to increase to gain energy from the 
field. We estimate the tunnelling rate out of these 
metastable states if the electric field E is weak. They 
are given by P ^ ^-Sgin exponential accuracy, where 
S*Q is the action of a bounce. This is the instanton so- 
lution that corresponds to the minimal action needed to 
go between the two minima and back [l^ |26| . Such an 
instanton has the shape of a bubble of typical size in 
the space direction and L^- in the time direction, and is 
schematically represented in Fig. If we denote «,j the 
coordinate in space and time respectively, then 



(8) 



where 



is the deviation from the ground state. We 



consider unit instantons with rriij — 1 inside the bubble, 
and rriij — outside. The region where (juj interpolates 
between these two values is the wall which encircles the 
bubble, which is very thin in the large A* = I limit con- 
sidered here. It is useful to recall that in the pure Sinc- 
Gordon model (obtained here taking all Q = 0, n'f = nio) 
the bounce instanton solution (with zero friction coeffi- 
cient) is a circle in the x, r plane, since the theory is 
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FIG. 1: An instanton starting from the ground state where 

1 (see text). 



Lorentz invariant. Here, the surface tension of the in- 
stanton waUs is highly anisotropic. For a "time-hke" waU 
paraUel to the x axis, the surface tension is the same as 
pure Sine-Gordon (with renormahzed parameters) since 
the disorder is time independant. The corresponding cost 
in the action is arL^ where the hne tension of such waUs 
is (Tt = 2tt/ {K*a). For a "space-hke" wall parallel to 
the r axis the surface tension a^ii) is a random vari- 
able. In first approximation the typical instanton now 
has a rectangular shape, bounded in x by two vertical 
segments parallel to the r-axis at coordinates x = io and 
X = io + Lx, chosen as places where crx{i) is small. The 
rectangle is closed by two "time-like" segments at r = tq 
and T = To + Lt-. 

Let us consider a segment of length Lr of the wall 
parallel to the r direction between sites i and i -I- 1. The 
extra action due to the presence of the instanton is: 
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Thus, for a unit instanton mi+i j = 1, the line tension 



TP = o'x(i)-^r: depends on space position i: 
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(10) 



with Qi — [fi] — fi + 1/2. One easily sees that gi is 
uniformly distributed on the interval [0, 1]. In particular 
there is a finite weight around gi = which corresponds 
to "weak points" in the construction of the ground state 
where one can bifurcate from point i up to the boundary 
at low energy cost to a state where the phase is shifted 
by 2tt on the right of i (or conversely — 27r on the left of 
i for the wall on the right). Although these states can 
be close in energy, the tunnelling rate to them is zero. 
To obtain a non zero tunnelling rate one must consider 



"a kink" i.e. tunnelling to a neighboring state where the 
phase is shifted by 27r between two walls. This is the 
tunnelling process described by the above instanton. 

The total action cost of the above rectangular instan- 
ton is thus: 
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A'KeLxLr 
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(11) 



Since the two smallest numbers in a set of L^/a random 
numbers are typically of order a/L^ one can estimate 
(^x{io) + <yx[io + Lx) ~ One then easily estimate 

the line tension (|10|l and by minimizing the action 1)11(1 
get the optimal size for the instanton (for small e) 



This yields a decay rate: 

p ^ K'aV^ = exp[ 



L°P' = l/(2e) 
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(13) 



where we have introduced a characteristic energy scale 
A ~ ?i'u*/^ioc associated with the localization length. 
Note that m*/^ioc is the pinning frequency [ij. For a 
simple sine-Gordon theory the dependence is e"^*^/*^^^-* 
and Am is the Mott gap. This expression corresponds to 
Zcner tunnelling across the gap. 

Although the above analysis is expected to give cor- 
rectly the electric field dependence, the precise prefac- 
tor in the exponential might be modified by additional 
physical effects and its precise determination, beyond 
the crude estimate given here, is delicate. First strictly 
speaking, in order to reach a stationary state some 
amount of dissipation should be introduced in the model. 
This dissipation changes the cost of the time variation of 
the phase and thus Ut but does not affect (Jx(i)- It thus 
slightly changes the prefactor which could in principle be 
studied as in Next since Uxiio) ^ cFx{io + Lx), the 
instanton has a lozenge shape and the space like portion 
can improve its action by taking advantage locally of fa- 
vorable pins. That may slightly renormalize downwards 
a-T- Let us also point out that to obtain the response 
of the system we have computed here a typical instan- 
ton, which can occur repeatedly in the volume of the 
system. There are rarer events that correspond to faster 
tunnelling. Let us divide the system in intervals of scale 
Rx ^ (up to the system size) . Within each interval 
there is typically one place to put two walls separated by 
L°P* and for which {(Jx{io) + (Tx{io + L°J>^) ~ Thus 
the ground state tunnels (back and forth) with these 
states at a much faster rate. However since these tun- 
nelling events correspond to special places the density of 
such atypical kinks being 0{L'^^ / Rx) they cannot lead to 
a macroscopic current. The system thus stays essentially 
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FIG. 2: Left: a bounce for large field (dashed line) and at very 
small fields (solid line). Since the bounce reaches the finite 
size in imaginary time due to the finite temperature, it opens 
up. Right : The growth of the action at zero temperatures 
(dashed line) has a maximum which depends on the electric 
field. This is the barrier to pass to do the tunnelling. At fi- 
nite temperature, the barrier decreases very rapidly when the 
bounce opens up (solid line). The barrier is thus determined 
essentially by the temperature. 



blocked until the tunnelling events due to the typical in- 
stantons can take place. It would however be interesting 
to see whether such rare events could serve as nucleation 
center for "quantum avalanches" , which could only in- 
crease the creep rate. Although this clearly goes beyond 
the present study, the explicit construction of the low ly- 
ing states presented here may allow for a precise study 
of this faster dynamics, at least numerically. 

Let us now see how this quantum creep which corre- 
sponds to tunnelling due to quantum effects is modified 
by the presence of a non zero temperature. Because of 
the finite temperature the time integration in imaginary 
time is limited to the finite value Lr,M = u*/3?i/(^ioc/a) 
(because of the rescaling). This means that the above 
analysis which was done at /3 = cx) remains valid as long 
as the size of the bounce in the time direction is smaller 
than ir,j\/- When the size of the bounce reaches the 
boundary the instanton opens (there are periodic bound- 
ary conditions in imaginary time) as shown in Fig. |21 
Because there is now no contribution coming from in- 
stantons parallel to the space direction, it is easy to see 
from that the action decreases linearly with the size 
Lx of the instanton. The tunnelling rate is thus fixed by 
the maximum barrier, i.e. the value of the action when 
the bounce reaches Lr,M- Because now the maximum 
barrier is not fixed by the electric field any more, one has 
to consider both the forward and backward jumps as for 
the standard TAPE argument j^. The net probability 
current is thus proportional to 
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(14) 



die point ifT^ when = Lr.M- One thus recovers below 
the crossover field e = ^ioc/(2/3?iau*) a linear response, 
with a conductivity proportional to 



(j{T) cx e 



exp[ 



(15) 



Quite remarkably the temperature dependence of the 
conductivity as obtained by the present formula is identi- 
cal to Mott's variable range hopping where the tran- 
sition between localized states close in energy is provided 
by external source of inelastic scattering such as the elec- 
tron phonon interaction. The important difference be- 
tween our result and the standard VRH law is that here, 
inelastic processes are coming from the electron-electron 
interaction itself (hidden in the existence of the Luttinger 
liquid parameter K) . Thus the pref actor in the exponen- 
tial contains electronic energy scales. The VRH formula 
contains normally the Debye temperature for phonons. 
Our result thus lead to a quite different energy scale in 
the exponential. Although our calculation is done in one 
dimension only, it is most likely that in higher dimen- 
sion as well one can obtain similar formulas. Let us note 
that in one dimension the above instanton picture is very 
similar physically to the VRH picture, if one remembers 
that in a Luttinger liquid a kink in (j) is related to the 
presence of a charge though the formula p = — V^/tt. 
Shifting the ground state by one unit is equivalent to 
moving an electron. 

TG would like to thank B. L. Altshuler for interesting 
discussions. 



where S* is the action of the bounce as given by the sad- 
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